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I. INTRODUCTION 



The first observation of meson by the CDF collaboration [|IJ] opens a new direction in the 
physics of hadrons containing heavy quarks. This particle completes the list of heavy quarkonia 
(QQ') and heavy flavored mesons accessible for the experimental investigations. It begins another 
list of long-lived hadrons containing two heavy quarks. So, in addition to the B+ meson this class 
of hadrons would be continued by the doubly heavy baryons S cc , E bc and E bb (see notations in 
the framework of quark model by PDG 0). The experimental discovery of was prepared by 
theoretical studies of the meson spectroscopy as well as the mechanisms of its production and decay 
(see review in 0). To observe the doubly heavy baryons it is necessary to give reliable theoretical 
predictions of their properties. The initial steps forward such the goal were done. 

1. In ref. Q the estimates for the lifetimes of 5+ and S+ + baryons were obtained in the framework 
of Operator Product Expansion in the inverse heavy quark mass. 

2. Papers || were devoted to the investigation of differential and total cross-sections for the 
production of Hqq/ baryons in various interactions in the model of fragmentation, in the model 
of intrinsic charm || (for the hadronic production of S cc ) and in the framework of perturbative 
QCD calculations up to 0(af) contributions, taking into account the hard nonfragmentational 
regime in addition to the fragmentation, which dominates at high transverse momenta pr ^> 
M. 

3. In refs. |7|,|8| the masses of ground states of doubly heavy baryons were estimated, and the 
excitations of E cc were considered in ref. 

In the present paper we analyze the basic spectroscopic characteristics for the families of doubly 
heavy baryons E QQ/ = (QQ'q), where q = u, d and VLqqi = (QQ's). 

A general approach of potential models to calculate the masses of baryons containing two heavy 
quarks was considered in refs. 0. The physical motivation used the pair interactions between the 
quarks composing the baryon, that was explored in the three-body problem. Clear implications for 
the mass spectra of doubly heavy baryons were derived. So, for the given masses of heavy charmed 
and beauty quarks, the approximation of factorization for the motion of doubly heavy diquark and 
light quark is not accurate. It results in the ground state mass and excitation levels, essentially 
deviating from the estimates in the framework of appropriate three-body problem. For example, we 
can easily find that in the oscillator potential of pair interactions an evident introduction of Jacobi 
variables leads to the change of vibration energy uj — > J\oj in comparison with naive expectations 
of diquark factorization. 

However, the string-like picture of doubly heavy baryon shown in Fig.|T], certainly destroys the 
above conclusions based on the pair interactions. Indeed, to the moment we have to introduce the 
centre of string, which is very close to the centre of mass for the doubly heavy diquark. Furthermore, 
the light quark interacts with the doubly heavy diquark as a whole, i.e. with the string tension 
identical to that in the heavy-light mesons. Therefore, two different assumptions on the nature 
of interactions inside the doubly heavy baryons: pair interactions or string-like picture, result 
in a distinct variation of predictions on the mass spectra of these baryons: the ground states 
and excitation levels. The only criterion testing the assumptions is provided by an experimental 
observation and measurements. 
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In this paper we follow the approximation of double heavy diquark, which is quite reasonable as 
we have clarified in the discussion given above. To enforce this point we refer to the consideration 
of doubly heavy baryon masses in the framework of QCD sum rules [[|, that recently was essentially 
improved in paper flQ| , exploring the NRQCD version of sum rule method. The conclusion drawn 



in the sum rules is the ground state mass is in a good agreement with the estimates obtained in the 
potential approach with the factorization of doubly heavy diquark. 

The qualitative picture for the forming of bound states in the system of {QQ'q) is determined 
by the presence of two scales of distances, which are given by the size of QQ'-diquark subsystem, 
tqqi, in the anti-triplet color state as well as by the confinement scale, Aq C d, for the light quark q, 
so that 

rQQi ■ Aqcd < 1, A QCD < m Q . 

Under such conditions, the compact diquark QQ' looks like a static source approximated by the 
local colored QCD field interacting with the light quark. Therefore, we may use a set of reliable 
results in models of mesons with a single heavy quark, i.e. with a local static source belonging 
to the anti-triplet representation of SU(3) C group. The successful approaches are the potential 
models |E| and the Heavy Quark Effective Theory (HQET) |13j in the framework of expansion 
in the inverse heavy quark mass. We apply the nonrelativistic quark model with the potential 
by Buchmiiller-Tye |TJ[ . Then theoretically we can talk on the rough approximation for the light 
quark. Indeed, since m® CD <C Aqcd the light quark is relativistic. Nevertheless, we introduce the 
system with a finite number of degrees of freedom and an instantaneous interaction V(r). This fact 
is a disadvantage because the confinement supposes the following: a) the generation of sea around 
the light quark, i.e. the presence of infinite number of gluons and quark-antiquark pairs, and b) 
the nonperturbative effects with the correlation time tqcd ~ 1/Aqod> that is beyond the potential 
approach. However, phenomenologically the introduction of constituent mass m^ p ~ Aqcd as a 
basic parameter determining the interaction with the QCD condensates, allows us successfully to 
adjust the nonrelativistic potential model with a high accuracy (8M ^ 30 — 40 MeV) by fitting the 
existing experimental data, that makes the approach to be quite a reliable tool for the prediction 
of masses for the hadrons, containing the heavy and light quarks. 

As for the diquark QQ', it is completely analogous to the heavy quarkonium QQ' except the 
very essential peculiarities. 

1. {QQ')z c is a system with the nonzero color charge. 

2. For the quarks of the same flavor Q = Q' it is necessary to take into account the Pauli principle 
for the identical fermions. 

The second item turns out to forbid the sum of quark spins S=0 for the symmetric, spatial parity P- 
even wave functions of diquark, ^(r) (the orbital momentum equals = 2n, where n — 0, 1, 2 . . . 
), as well as S=l is forbidden for the anti-symmetric, -odd functions ^d(r) (i.e. = In + 1). The 
nonzero color charge leads to two problems. 

First, we cannot generally apply the confinement hypothesis on the form of potential (an infinite 
growth of energy with the increase of the system size) for the object under consideration. However, 
it is unpossible to imagine a situation, when a big colored object with a size r > 1/Aqcd has a finite 
energy of self-action, and, to the same moment, it is confined inside a white hadron (the singlet over 
SU(3) C ) with r ~ 1 /Aqcd due to the interaction with another colored source. In the framework 
of well-justified picture of the hadronic string, the tension of such string in the diquark with the 
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external leg inside the baryons is only two times less than in the quark-antiquark pair inside the 
meson qq', and, hence, the energy of diquark linearly grows with the increase of its size. So, the effect 
analogous to the confinement of quarks takes place in the similar way. In the potential models we 
can suppose that the quark binding appears due to the effective single exchange by a colored object 
in the adjoint representation of SU(3) C (the sum of scalar and vector exchanges is usually taken). 
Then, the potentials in the singlet (qq') and anti-triplet (qq') states differ by the factor of 1/2, that 
means the confining potential with the linear term in the QCD-motivated models for the heavy 
diquark (QQ')%- In the present work we use the nonrelativistic model with the Buchmiileer-Tye 
potential for the diquark, too. 

Second, in the singlet color state (QQ') there are separate conservations of the summed spin 
S and the orbital momentum L, since the QCD operators for the transitions between the levels 
determined by these quantum numbers, are suppressed. Indeed, in the framework of multipole 



expansion in QCD [15[], the amplitudes of chromo-magnetic and chromo-electric dipole transitions 
are suppressed by the inverse heavy quark mass, but in addition, the major reason is provided by 
the following: a) the necessity to emit a white object, i.e. at least two gluons, which results in the 
higher order in 1/tjiq, and b) the projection to a real phase space in a physical spectrum of massive 
hadrons in contrast to the case of massless gluon. Furthermore, the probability of a hybrid state, 
say, the octet sybsystem (QQ') and the additional gluon, i.e. the Fock state \QQ' 8c g), is suppressed 
due to both the small size of system and the nonrelativistic motion of quarks (for a more strict 



consideration see ref. [jrjj). In the anti-triplet color state, the emission of a soft nonperturbative 



gluon between the levels determined by the spin Sd and the orbital momentum Ld of diquark, is 
not forbidden, if there are no some other no-go rules or small order-parameters. For the quarks 
of identical flavors inside the diquark, the Pauli principle leads to that the transitions are possible 
only between the levels, which either differ by the spin (ASd = 1) and the orbital momentum 
(ALd = 2n + 1), instantaniously, or belong to the same set of radial and orbital excitations with 
ALd = In. Therefore, the transition amplitudes are suppressed by a small recoil momentum of 
diquark in comparison with its mass. The transition operator changing the diquark spin as well as 
its orbital momentum, has the higher order of smallness because of either the additional factor of 
1/m.Q or the small size of diquark. These suppressions lead to the existence of quazi-stable states 
with the quantum numbers of Sd and Ld- In the diquark composed by the quarks of different flavors, 
6c, the QCD operators of dipole transitions with the single emission of soft gluon are not forbidden, 
so that the lifetimes of levels can be comparable with the times for the forming of bound states or 
with the inverse distances between the levels themselves. Then, we cannot insist on the appearance 
of excitation system for such the diquark with definite quantum numbers of the spin and orbital 
momentum^. 

Thus, in the present work we explore the presence of two physical scales in the form of factoriza- 
tion for the wave functions of the heavy diquark and light constituent quark. So, in the framework 
of nonrelativistic quark model the problem on the calculation of mass spectrum and characteristics 
of bound states in the system of doubly heavy baryon is reduced to two standard problems on the 
study of stationary levels of energy in the system of two bodies. After that, we take into account the 



*In other words, the presence of gluon field inside the baryon Efe c leads to the transitions between the 
states with the different excitations of diquark, like \bc) — ► \bcg) with ASd = 1 or ALd = lj which are not 
suppressed. 
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relativistic corrections dependent of the quark spins in two subsystems under consideration. The 
natural boundary for the region of stable states in the doubly heavy system can be assigned to the 
threshold energy for the decay into a heavy baryon and a heavy meson. As was shown in ||17||, the 



appearance of such threshold in different systems can be provided by the existence of an universal 
characteristics in QCD, a critical distance between the quarks. At distances greater than the critical 
separation, the quark-gluon fields become unstable, i.e. the generation of valence quark-antiquark 
pairs from the sea takes place. In other words, the hadronic string having a length greater than 
the critical one, decay into the strings of smaller sizes with a high probability close to unit. In the 
framework of potential approach this effect can be taken into account by that we will restrict the 
consideration of excited diquark levels by the region, wherein the size of diquark is less than the 
critical distance, tqqi < r c « 1.4 — 1.5 fm. Furthermore, the model with the isolated structure of 
diquark looks to be reliable, just if the size of diquark is less than the distance to the light quark 
r QQ' < n- 

The peculiarity of quark-diquark picture for the doubly heavy baryon is the possibility of mixing 
between the states of higher diquark excitations, possessing the different quantum numbers, because 
of the interaction with the light quark. Then it is difficult to assign some definite quantum numbers 
to such excitations. We will discuss the mechanism of this effect. 

The paper is organized as follows. In Section II we describe a general procedure for the calcu- 
lation of masses for the doubly heavy baryons in the framework of assumptions drawn above. We 
take into account the spin-dependent corrections to the potential motivated in QCD. The results of 
numerical estimates are presented in Section III, and, finally, our conclusions are given in Section 
IV. 



II. NONRELATIVISTIC POTENTIAL MODEL 

As we have mentioned in the Introduction, we solve the problem on the calculation of mass 
spectra of baryons containing two heavy quarks, in two steps. First, we compute the energy levels 
of diquark. Second, we consider the two body problem for the light quark interacting with the 
point-like diquark having the mass obtained in the first step. In accordance with the effective 
expansion of QCD in the inverse heavy quark mass, we separate two stages of such the calculations. 
So, the nonrelativistic Schrodinger equation with the model potential motivated by QCD, is solved 
numerically. After that, the spin- dependent corrections are introduced as perturbations suppressed 
by the quark masses. 



A. Potential 

We use the Buchmiiller-Tye potential, which takes into account the asymptotic freedon of QCD 
at short distances. So, the effective coupling constant in the exchange by the octet color state 
between the quarks is approximated by the QCD running coupling constant up to two-loop accuracy. 
At long distances, there is the linear term of interaction energy, which leads to the confinement. 
These two regimes are the limits for the effective /3-function by Gell-Mann-Low. It was given 



explicitly in [14]. In the anti-triplet quark state we introduce the factor of 1/2 because of the color 
structure of bound quark-quark state. For the interaction of diquark with the light constituent 
quark, the corresponding factor is equal to unit. 
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As was shown in ||18|| , the nonperturbative constituent term introduced into the mass of nonrela- 
tivistic quark, exactly coincides with the additive constant, subtracted from the coulomb potential. 
We extract the masses of heavy quarks by fitting the real spectra of charmonium and bottomonium, 



m c = 1.486 GeV, m b = 4.88 GeV, (1) 

so that the mass of the level in the heavy quarkonium has been calculated as, say, M(c) = 2m c + E, 
where E is the energy of stationary Schrodinger equation with the model potential V. Then, we have 
supposed that the mass of meson with a single heavy quark is equal to M(Qq) = tuq + m q + E, and 
E = (T) + (V — SV), whereas the additive term in the potential is introduced because the constituent 
mass of light quark is determined as a part of interaction energy 5V = m q . In accordance with 
fitting the masses of heavy mesons, we get m q = 0.385 GeV. 

The results of calculations for the eigen-energies in the Schrodinger equation with the 
Buchmiiller-Tye potential are presented in Tabs. 1-3. The characteristics of corresponding wave 
functions are given in Tabs. 4-6. 

We have checked that with a good accuracy the binding energy and the wave function of light 
quark practically do not depend on the flavors of heavy quarks. Indeed, large values of diquark 
masses give small contributions into the reduced masses. This fact leads to small corrections to 
the wave functions in the Schrodinger equation. So, for the states lying below the threshold of 
doubly heavy baryon decay into the heavy baryon and heavy meson, the energies of levels of light 
constituent quark are equal to 

E(ls) = 0.38 GeV, E{2s) = 1.09 GeV, E(2p) = 0.83 GeV, 

where the energy has been defined as the sum of light quark constituent mass and eigen-value of 
Schrodinger equation. In HQET the value of A = E(ls) is generally introduced. Then we can draw 
a conclusion that our estimate of A is in a good agreement with calculations in other approaches 
|i~3| . This fact confirms the reliability of such the phenomenological predictions. For the light quark 
radial wave functions at the origin we find 

R u (0) = 0.527 GeV 3/2 , R 2s (0) = 0.278 GeV 3/2 , R' 2p (0) = 0.127 GeV 5/2 . 

The analogous characteristics of bound states of the c-quark interacting with the 66-diquark, are 
equal to 

E(ls) = 1.42 GeV, E(2s) = 1.99 GeV, E(2p) = 1.84 GeV, 
with the wave functions 

R l8 (0) = 1.41 GeV 3/2 , R 2s (0) = 1.07 GeV 3/2 , R' 2p (0) = 0.511 GeV 5/2 . 

For the binding energy of strange constituent quark we add the current mass m s « 100 — 150 MeV. 



B. Spin-dependent corrections 



According to [13] , we introduce the spin-dependent corrections causing the splitting of nL-levels 
of diquark as well as in the system of light constituent quark and diquark (n = n r + L + 1 is the 
principal number, n r is the number of radial excitation, L is the orbital momentum). For the heavy 
diquark containing the identical quarks we have 
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{d) 1 /Ld_Sd\ / dV(r) 8 1 \ 

~\ a *\ AT 2 o t 6 ( L d • S d) 2 + 3(L d • S d ) - 2L d 2 S d 2 ]^, 
3 4L d - 3 r 3 

where L d , S d are the orbital momentum in the diquark system and the summed spin of quarks 
composing the diquark, respectively. Taking into account the interaction with the light constituent 
quark gives (S = S d + Si) 

(0 1 /L ■ S d 2L-SA f_dV{r) 8 1 \ 

^w-il^l + 2m 2 )y rdr +3 a v3 y ) 

+k 1 ( L - S <* + 2L - S + jk-J— (S d + L d ) ■ Sl [47rf(r)] (3) 
3 rnqrui r 6 3 omgrrii 

1 AT2 X , [6(L ■ S) 2 + 3(L ■ S) - 2L 2 S 2 
3 tjiqitii 4Lr — 3 

-6(L-S d ) 2 -3(L-S d ) + 2L 2 S d 2 ]^, 

where the first term corresponds to the relativistic correction to the effective scalar exchange, and 
other terms appear because of corrections to the effective single-gluon exchange with the coupling 
constant a s . 

The value of effective parameter a s can be determined in the following way. The splitting in the 
S"-wave heavy quarkonium (Q1Q2) is given by the expression 

AM(ns) = la s —^—\R nS (0)\ 2 , (4) 
9 mim 2 

where R n s{r) is the radial wave function of quarkonium. From the experimental data on the system 
of cc 

AM(1S, cc) = 117 ± 2 MeV, (5) 

and R\s(0) calculated in the model, we can determine a s (ty). Let us take into account the de- 
pendence of this parameter on the reduced mass of the system, fi. In the framework of one-loop 
approximation for the running coupling constant of QCD we have 

whereas 6 = 11 — 2nf/3 and rif = 3 at p 2 < m 2 . From the phenomenology of potential models [0 
we know that the average kinetic energy of quarks in the bound state practically does not depend 
on the flavors of quarks, and it is given by the values 

(T d ) « 0.2 GeV, (7) 
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and 



(T t ) » 0.4 GeV, (8) 

for the anti-triplet and singlet color states, correspondingly. Substituting the definition of the 
nonrelativistic kinetic energy 

(T) = |>, (9) 

we get 

Air 

a - (p2) = fr.ln(2( WCT r <10) 

whereas numerically Aqcd ~ 113 MeV. 

For the identical quarks inside the diquark, the scheme of LS'-coupling well known for the 
corrections in the heavy quarkonium, is applicable. Otherwise, for the interaction with the light 
quark we use the scheme of jj-coupling (here, LSi is diagonal at the given Ji, (Ji = L+Si, J = Ji+J), 
where J denotes the total spin of baryon, and J is the total spin of diquark, J = Sd + Ld)- 

Then, to estimate various terms and mixings of states, we use the transformations of bases (in 
what follows S = Si + J) 



\J; Ji) = E(-i) (J+5i+L+J) y/(?s + 1)(2 j, + 1) 1 1 S J 1}\J;S), (ii) 



and 



J; Ji) = E(-i) (J+Si+L+J V( 2 ^ + i)( 2 ^ + i){ s t j 3 j\ } 1 J; Jd 



(12) 



Thus, we have defined the procedure of calculations for the mass spectra of doubly heavy baryons. 
This procedure leads to results presented in the next section. 



III. NUMERICAL RESULTS 



In this Section we present the results on the mass spectra with account for the spin-dependent 
splitting of levels. As we have clarified in the Introduction, the doubly heavy baryons with identical 
heavy quarks allow quite a reliable interpretation in terms of diquark quantum numbers (the summed 
spin and the orbital momentum). Dealing with the excitations of frc-diquark, we show the results 
on the spin-dependent splitting of the ground lS-state, since the emission of soft gluon breaks the 
simple classification of levels for the higher excitations of such diquark. 

For the doubly heavy baryons, the quark-diquark model of bound states obviously leads to the 
most reliable results for the system with the larger mass of heavy quark, i.e. for E bb . 
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A. baryons 

For the quantum numbers of levels, we use the notations UdLdfiili, i.e. we show the value of 
principal quantum number of diquark, its orbital momentum by a capital letter and the principal 
quantum number for the excitations of light quark and its orbital momentum by a lower-case letter. 
We denote the shift of level by A( J ) as dependent on the total spin of baryon J. So, for \S2p we 
have 

A ( §) = 10.3 MeV. (13) 

The states with the total spin J = | (or |), can have different values of Ji, and, hence, they have a 
nonzero mixing, when we perform the calculations in the perturbation theory built over the states 
with the definite total momentum J\ of the light constituent quark. For J = |, the mixing matrix 
equals 

so that the mixing practically can be neglected, and the level shifts are determined by the values 

A /( i } = A; = -3.0 MeV, (15) 
A ( i) = Ai = 11.4 MeV. 



For J = |, the mixing matrix has the form 



-5.7 -17.8 
-17.8 -14.9 



MeV, (16) 



with the eigen-vectors given by 



|lS2p(l')) = 0.790| J, = \) - 0.613| J, = I), (17) 
\\S2p(\)) = 0.613| J ; = -) + 0.790| J t ' 



2' 



and the eigen-values equal 



A' ( ^) = A' 2 = 8.1 MeV, (18) 
A ( 5) = A 2 = -28.7 MeV. 

For the 2S'2p-level, the corresponding quantities are equal to 

A ( f) = 10.3 MeV, (19) 

and for J = |, the mixing matrix is equal to 
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so that 

A'(§) = A; = -3.6 MeV, (21) 
A ( i) = Ax = 12.4 MeV. 

For J = |, the matrix has the form 

/ -6.1 -17.6' 



-17.6 -13.5 I MeV ' < 22 > 



with the eigen-vectors 



\lS2p(^)) = 0.776| J, = ^) - 0.631| J, = ^>, (23) 

|152p(i)) = 0.631| J, = |) + 0.776| J z = ±), 
possessing the eigen-values 

A' ( ^) = A' 2 = 8.2 MeV, (24) 
A ( ^ = A 2 = -27.8 MeV. 

We can straightforwardly see, that the difference between the wave functions as caused by the 
different masses of diquark subsystem, is unessential. 
The splitting of diquark, A^ Jd \ has the form 

3Dls : 
A (3) = -0.06 MeV, 

A (2) = 0.2 MeV, (25) 
A (1) = -0.2 MeV. 

4£>ls : 
A (3) = -2.6 MeV, 

A (2) = -0.8 MeV, (26) 
A (1) = -4.6 MeV. 

5Dls : 
A (3) = 2.6 MeV, 

A (2) = -0.9 MeV, (27) 
A (1) = -4.7 MeV. 

5Gls : 
A (5) = -0.3 MeV, 
A {4) = 0.3 MeV, 

A (3) = 1.1 MeV, (28) 
A (2) = 1.7 MeV, 
A (1) = 2.0 MeV. 
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QGls : 
A (5) = 3.2 MeV, 
A (4) = -0.5 MeV, 

A (3) = -4.4 MeV, (29) 
A< 2 > = -7.9 MeV, 
A (1) = -10.5 MeV. 

Such the corrections are unessential up to the current accuracy of method (SM ^ 30 — 40 MeV). 
They can be neglected for the excitations, whose sizes are less than the distance to the light quark, 
i.e. for the diquarks with small values of principal quantum number. 

For the hyper-fine spin-spin splitting in the system of quark-diquark, we have 



A?\ = - 



J(J+1)_ J(J+1)-| 



a s (2 f iT)^—\R l (0)\ 2 , (30) 
m c mi 



where Ri(0) is the radial wave function at the origin for the light constituent quark, and for the 
analogous shift of diquark level, we find 

Ag>. = l -a s {2^T)^-\RM\ 2 - (31) 

The mass spectra of and baryons are shown in Fig.l and in Tab. 7, wherein we restrict 
ourselves by the presentation of S-, P- and D-wave levels. 

We can see in Fig.l, that the most reliable predictions are the masses of baryons 151s (J p = 
3/2+, 1/2+), 2Pls (J p = 3/2", 1/2") and 3Dls {J p = 7/2+, . . . 1/2+). The 2Pls-level is quazi- 
stable, because the transition into the ground state requires the instantaneous change of both 
the orbital momentum and the summed spin of quarks inside the diquark. The analogous kind 
of transitions seems to be the transition between the states of ortho- and para-hydrogen in the 
molecule of H 2 . This transition take place in a non- homogeneous external field due to the magnetic 
moments of other molecules. For the transition of 2Pls — > 151s, the role of such the external field 
is played by the non-homogeneous chromo-magnetic field of the light quark. The corresponding 
perturbation has the form 

5V ~ —[Si ■ Hx + S 2 • H 2 - (Si + S 2 ) • (H)] 

m Q 

= J-(V • r d ) (Si - S 2 ) • H ~ — (Si - S 2 ) • J, fin), 

2m Q m Q m q rf 

where f(ri) is a dimensionless nonperturbative function depending on the distance between of 
the light quark and diquark. The 5V operator changes the orbital momentum of light quark, 
too. It results in the mixing between the states with the same values of J p . If the splitting 
is not small (for instance, 2Pls — \S2p, where AE ~ Aq C d), then the mixing is suppressed as 
5V/AE ~ -^tJ^^e ^ L Since tne admixture of \S2p in the 2Pls-state is low, the 2Pls-levels 
are quazi-stable, i.e. their hadronic transitions into the ground state with the emission of 7r-mesons 
are suppressed as we have derived, though an additional suppression is given by a small value of 
phase space. Therefore, we have to expect the presence of narrow resonances in the mass spectra 
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of pairs H&&7T, as they are produced in the decays of quazi-stable states with J p = 3/2~~, 1/2 - . The 
experimental observation of such levels could straightforwardly confirm the existence of diquark 
excitations and provide the information on the character of dependency in f(ri), i.e. on the non- 
homogeneous chromo-magnetic field in the nonperturbative region. 

Sure, the 3_Dls J p = 7/2 + , 5/2 + states are also quazi-stable, since in the framework of multipole 
expansion in QCD they transform into the ground state due to the quadru-pole emission of gluon 
(the E2-transition with the hadronization gq — > q'n) . 

As for the higher excitations, the 3Pls-states are close to the lS2p-levels with J p = 3/2 _ , l/2 _ , 
so that the operators changing both the orbital momentum of diquark and its spin, can lead to the 
essential mixing with an amplitude SV nn r/AE nn i ~ 1, despite of suppression by the inverse heavy 
quark mass and small size of diquark. We are sure that the mixing slightly shifts the masses of 
states. The most important effect is a large admixture of lS2p in 3Pls. It makes the state to 
be unstable because of the transition into the ground 1 Sis-state with the emission of gluon (the 
El-transition). This transition leads to decays with the emission of 7r-mesons[]. 

The level lS2p J p = 5/2 _ has the definite quantum numbers of diquark and light quark motion, 
because there are no levels with the same values of J p in its vicinity. However, its width of transition 
into the ground state and 7r-meson is not suppressed and seems to be large, T ~ 100 MeV. 

The following transitions take place 
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The D-wave transitions are suppressed by the ratio of low recoil momentum to the mass of baryon. 

The width of state J p = 3/2 + is completely determined by the radiative electromagnetic Ml- 
transition into the basic J p = l/2 + state. 

B. H cc baryons 

The calculation procedure described above leads to the results for the doubly charmed baryons 
as presented below. 

For lS2p, the splitting is equal to 

A ( f ) = 17.4 MeV. (32) 
For J = |, the mixing is determined by the matrix 



^Remember, that the SQQ'-baryons are the iso-dublets. 
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4.3 -1.7 \ 
-1.7 7.8 MeV ' 



so that the eigen- vectors 



|lS2p(^)> = 0.986| J, = ^) + 0.164| J, = ±), 



3,, 3, 1 

2'' 



|lS2p(-)) = —0.164| J| = -) + 0.986| J, 



have the eigen- values 



A /( 2 ) = A; = 3.6 MeV, 
A ( i } = Ai = 8.5 MeV. 



For J —\i the mixing matrix equals 



where the vectors 



have the eigen- values 



-3.6 -55.0 
-55.0 -73.0 ) ' 



|lSr2p(|')> = 0.957| Ji = |> - 0.291| J, = I), 
\\S2 V <X)) = 0.291| J, = h + 0.957| J, = 



A ,( 5) = A' 2 = 26.8 MeV, 
A ( ^) = A 2 = -103.3 MeV. 

The splitting of 3D diquark level is given by 

A {3) = -3.02 MeV, 
A {2) = 2.19 MeV, 
A (1) = 3.39 MeV. 

Further, we have to take into account the hyper-fine spin-spin corrections in the quark 
system. 

For the IS- and 2S'-wave levels of diquark, the shifts of vector states are equal to 

A(1S) = 6.3 MeV, 
A(2S) = 4.6 MeV. 

The mass spectra of the S+ + and S+ baryons are presented in Fig. 2 and Tab. 8. 
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C. E bc baryons 

As we have already mentioned in the Introduction, the heavy diquark composed of the quarks 
of different flavors, turns out to be unstable under the emission of soft gluons. So, in the Fock state 
of doubly heavy baryon, there is a sizable nonperturbative admixture of configurations including 
the gluons and diquark with the various values of its spin and orbital momentum L d 

\B bcq ) = B \bc S 3 d c ' Ld , q) + H^bclf 1 '^, g, q) + H 2 \bc S £ Ld± \ g, <?> + ..., 

whereas the amplitudes of Hi, H 2 are not suppressed with respect to Ob. In the heavy quarkonium, 
the analogous operators for the octet-color states are suppressed by the probability of emission by 
the nonrelativistic quarks inside a small volume determined by the size of singlet-color system of 
heavy quark and antiquark. In the baryonic system under consideration, a soft gluon is restricted 
only by the ordinary scale of confinement, and, hence, there is no suppression. 

We suppose that the calculations of masses for the excited S bc baryons are not so justified in 
the given scheme. Therefore, we present only the result for the ground state with J p = l/2 + 

M*> = 6.85 GeV, M*. = 6.82 GeV, 

whereas for the vector diquark we have assumed that the spin-dependent splitting due to the 
interaction with the light quark is determined by the standard contact coupling of magnetic moments 
for the point-like systems. The picture for the baryon levels is shown in Fig. 3 with no account for 
the spin-dependent perturbations suppressed by the heavy quark masses. 

D. The doubly heavy baryons with the strangeness £Iqq*. 

In the leading approximation, we suppose that the wave functions and the excitation energies 
of strange quark in the field of doubly heavy diquark repeat the characteristics for the analogous 
baryons containing the ordinary quarks u, d. Therefore, the level system of baryons Qqq' reproduces 
that of Sqq/ up to an additive shift of the masses by the value of current mass of strange quark, 
m s « M(D S ) - M(D) « M(B S ) - M(B) « 0.1 GeV. 

Further, we suppose that the spin-spin splitting of 2Pls and 3-Dls levels of VIqqi is 20-30% 
less than in Sqq/ (the factor of m u ^/m s ). As for the 15'2p-level, the procedure described above 
can be applied. So, for Q bb , the matrix of mixing for the states with the different values of total 
momentum J\ practically can be assigned to be diagonal. This fact means that the following term 
of perturbation is dominant: 

1 / 2L- SA / dVjr) 8 1_\ 
4 y 2mf J \ rdr 3 V 3 y 

Therefore, we can think that the splitting of lS2p is determined by the factor of m 2 ud jvn 1 s with 
respect to the splitting of corresponding E bb , i.e. it is 40 % less than in E bb . Hence, the splitting is 
very small. 

For the baryon Q, the factor of m s /m c is not small. Hence, for \S2p, the mixing matrix is not 
diagonal, so that the arrangement of lS2p states of Q can be slightly different from that of S. 

The following peculiarity of VLqqi is of great interest: the low-lying S- and P-excitations of 
diquark are stable. Indeed, even after taking into account the mixing of levels, a gluon emission 
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makes a hadronization into the i^-meson (the transitions of Qqq' — > Sqq/ + K), while a single 
emission of 7r-meson is forbidden because of the conservation of iso-spin and strangeness. The 
hadronic transitions with kaons are forbidden because of insufficient splitting between the masses 
of VLqqi and Hqq/. The decays with the emission of pion pairs belonging to the iso-singlet state, are 
suppressed by a small phase space or even forbidden. Thus, the radiative electromagnetic transitions 
into the ground state are the dominant modes of decays for the low-lying excitations of VLqqi. 

E. Q bbc baryons 

In the framework of quark-diquark picture, we can build the model for the baryons containing 
three heavy quarks, bbc. However, as we estimate, the size of diquark turns out to be comparable 
with the average distance to the charmed quark. So, the model assumption on the compact heavy 
diquark cannot be quite accurate for the calculations of mass levels in this case. The spin-dependent 
forces are negligibly small inside the diquark, as we have already pointed out above. The spin-spin 
splitting of vector diquark interacting with the charmed quark, is given by 

A(ls) = 33 MeV, A(2s) = 18 MeV. 

For lS2p, the level shifts are small. So, for the state J p = l/2~ we have to add the correction 
of —33 MeV. For the 3-Dls-state the splitting is determined by the spin-spin interaction. The 
characteristics of excitations for the charmed quark in the model with the potential by Buchmiiller- 
Tye have been presented above. Finally, we obtain the picture of f2 bbc levels presented in Fig. 4 and 
Tab.9. 

Further, the excitations of ground Q bbc state can strongly mix with large amplitudes because of 
small splittings between the levels, but they have small shifts of masses. This effect takes place for 
3Pls - lS2p with J p = 1/2- , 3/2", and for 2Sls - 3Dls with J p = 1/2+ 3/2+. We suppose the 
prediction to be quite reliable for the states of 151s with J p = 1/2+, 3/2+, lS2p with J p = 5/2~ 
and 3-Dls with J p = 5/2+, 7/2+. For these excitations, we might definitely predict the widths 
of their radiative electromagnetic transitions into the ground state in the framework of multipole 
expansion in QCD. The widths for the transitions will be essentially determined by the amplitudes 
of admixtures, which have a strong model dependence. Therefore, the experimental study of elec- 
tromagnetic transitions in the family of Q bbc baryons could provide a valuable information on the 
mechanism of mixing between the different levels in the baryonic systems. The electromagnetic 
transitions combined with the emission of pion pairs, if not forbidden by the phase space, saturate 
the total widths of excited fl bbc levels. The characteristic value of total width is about T ~ 10 — 100 
keV, in the order of magnitude. 

Thus, the system of Q bbc can be characterized by a large number of narrow quazi-stable states. 

IV. CONCLUSION 

In this paper we have calculated the spectroscopic characteristics of baryons containing two heavy 
quarks, in the model with the quark-diquark factorization of wave functions. We have explored the 
nonrelativistic model of constituent quarks with the potential by Buchmuller-Tye. The region of 
applicability of such the approximations has been pointed out. 
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We have taken into account the spin-dependent relativistic corrections to the potential in the 
subsystems of diquark and light quark- diquark. Below the threshold of decay into the heavy baryon 
and heavy meson, we have found the system of excited bound states, which are quazi-stable under 
the hadronic transitions into the ground state. We have considered the physical reasons for the 
quazi-stability taking place for the baryons with two identical quarks. In accordance with the Pauli 
principle, the operators responsible for the hadronic decays and the mixing between the levels, are 
suppressed by the inverse heavy quark mass and the small size of diquark. This suppression is caused 
by the necessity of instantaneous change in both the spin and the orbital momentum of compact 
diquark. In the baryonic systems with two heavy quark and the strange quark, the quazi-stability 
of diquark excitations is provided by the absense of transitions with the emission of both a single 
kaon and a single pion. These transitions are forbidden because of small splitting between the levels 
and the conservation of the iso-spin and strangeness. 

The characteristics of wave functions can be used in calculations of cross-sections for the doubly 
heavy baryons in the framework of quark-diquark approximation. 

This work is in part supported by the Russian Foundation for Basic Research, grants 99-02-16558 
and 00-15-96645. The work of A.I.Onishchenko was supported, in part, by International Center of 
Fundamental Physics in Moscow, International Science Foundation, and INTAS-RFBR-95I1300 
grants. 
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1 ■ 11 1 

diquark level 


mass (GeV) 


/ 2\l/2 Ic \ 

(r z ) L/z (fm) 


1 ■ 11 1 

diquark level 


mass (GeV) 


; 9\l/2 IC \ 

\r ) / (fm) 


IS 


9.74 


0.33 


2P 


9.95 


0.54 


2S 


10.02 


0.69 


3P 


10.15 


0.86 


3S 


10.22 


1.06 


4P 


10.31 


1.14 


4S 


10.37 


1.26 


5P 


10.45 


1.39 


5S 


10.50 


1.50 


6P 


10.58 


1.61 


3D 


10.08 


0.72 


4D 


10.25 


1.01 


5D 


10.39 


1.28 


6D 


10.53 


1.51 


4F 


10.19 


0.87 


5F 


10.34 


1.15 


6F 


10.47 


1.40 


5G 


10.28 


1.01 


6G 


10.42 


1.28 


6M 


10.37 


1.15 



TABLE I. The spectrum of 66-diquark levels without spin-dependent splittings: masses and mean- 
squared radii. 



diquark level 


mass (GeV) 


(r 2 ) 1 / 2 (fm) 


diquark level 


mass (GeV) 


{r 2 ) 1/2 (fm) 


IS 


6.48 


0.48 


3P 


6.93 


1.16 


2S 


6.79 


0.95 


4P 


7.13 


1.51 


3S 


7.01 


1.33 


3D 


6.85 


0.96 


2P 


6.69 


0.74 


4D 


7.05 


1.35 


4F 


6.97 


1.16 


5F 


7.16 


1.52 


5G 


7.09 


1.34 


6H 


7.19 


1.50 



TABLE II. The spectrum of 6c-diquark levels without spin-dependent splittings: masses and mean- 
squared radii. 



diquark level 


mass (GeV) 


(r 2 ) 1 / 2 (fm) 


diquark level 


mass (GeV) 


(r 2 ) 1 / 2 (fm) 


IS 


3.16 


0.58 


3P 


3.66 


1.36 


2S 


3.50 


1.12 


4P 


3.90 


1.86 


3S 


3.76 


1.58 


3D 


3.56 


1.13 


2P 


3.39 


0.88 


4D 


3.80 


1.59 



TABLE III. The spectrum of cc-diquark levels without spin-dependent splittings: masses and 
mean-squared radii. 
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nL 




nL 




IS 


1.346 


2P 


0.479 


2S 


1.027 


3P 


0.539 


3S 


0.782 


4P 


0.585 


4S 


0.681 


5P 


0.343 



TABLE IV. The characteristics of radial wave function for the 66-diquark: R d ( ns )(0) (GeV 3 / 2 ), -R^( np )(0) 
(GeV 5 / 2 ). 



nL 




nL 


R'd(n P )(°) 


IS 




0.726 


2P 




0.202 


2S 




0.601 


3P 




0.240 


3S 




0.561 


4P 






TABLE V. The characteristics of radial wave function for the frc-diquark: R d ( ns 


(0) (GeV 3 / 2 ),^ (np) (0) 


(GeV 5 / 2 ). 








nL 


Rd(ns) (0) 


nL 




IS 




0.530 


2P 




0.128 


2S 




0.452 


3P 




0.158 



TABLE VI. The characteristics of radial wave function for the cc-diquark: R d ( n s){ty (GeV 3 / 2 ), R d ( np -)(0) 
(GeV 5 / 2 ). 



(n d L d mLi), J p 


mass (GeV) 


(n d L d niLi), J p 


mass (GeV) 


(IS ls)l/2+ 


10.093 


(3P Is) 1/2- 


10.493 


(IS ls)3/2+ 


10.133 


(3D ls)5/2'+ 


10.497 


(2P ls)l/2 _ 


10.310 


(3D ls)7/2+ 


10.510 


(2P ls)3/2 _ 


10.343 


(3P ls)3/2" 


10.533 


(2S ls)l/2+ 


10.373 


(1S 2p)l/2~ 


10.541 


(2S ls)3/2+ 


10.413 


(1S 2p)3/2~ 


10.567 


(3D ls)5/2+ 


10.416 


(1S 2p)l/2'- 


10.578 


(3D ls)3/2'+ 


10.430 


(1S 2p)5/2~ 


10.580 


(3D ls)l/2+ 


10.463 


(IS 2p)3/2'- 


10.581 


(3D ls)3/2+ 


10.483 


(3S ls)l/2+ 


10.563 



TABLE VII. The mass spectrum of E bb and E bb baryons. 
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/_, T „ T \ TP 

(n d L d niLi), J r 


mass (uevj 


{n d L d niLi), 


mass (GeV) 


(1b lsjl/z^ 


o.47o 


(3r Isjl/z 


3.97z 


(lb ls)6/Z^ 


3.61 


/or\ 1 \o /o/-(- 

(3D ls)3/2 T 


4.007 


/on -i \i /o — 

[Zr lsjl/2 


3.702 


(lb 2p)6/2 


a no A 

4.034 


/or 1 \r /0+ 

(3D Isj5/2 T 


3.781 


(lb 2p)3/2 


A non 

4.039 


f2S lsH/2+ 


3.812 


(IS 20)5/2" 


4.047 


(3D ls)3/2+ 


3.83 


(3D ls)5/2'+ 


4.05 


(2P ls)3/2- 


3.834 


(IS 2p)l/2'- 


4.052 


(3D ls)l/2+ 


3.875 


(3S ls)l/2+ 


4.072 


(IS 2p)l/2~ 


3.927 


(3D ls)7/2+ 


4.089 


(2S ls)3/2+ 


3.944 


(3P ls)3/2- 


4.104 


TABLE VIII. The mass spectrum of H+ + and baryons. 


(n d L d niLi), J p 


mass (GeV) 


(n d L d niLi), J p 


mass (GeV) 


(IS ls)l/2+ 


11.12 


(3D ls)3/2'+ 


11.52 


(IS ls)3/2+ 


11.18 


(3D ls)5/2'+ 


11.54 


(2P ls)l/2 _ 


11.33 


(IS 2p)l/2~ 


11.55 


(2P ls)3/2 - 


11.39 


(3D ls)7/2+ 


11.56 


(2S ls)l/2+ 


11.40 


(IS 2p)3/2'- 


11.58 


(3D ls)5/2+ 


11.42 


(1S 2p)3/2~ 


11.58 


(3D ls)3/2+ 


11.44 


(IS 2p)l/2'- 


11.59 


(3D ls)l/2+ 


11.46 


(1S 2p)5/2- 


11.59 


(2S ls)3/2+ 


11.46 


(3P ls)3/2 _ 


11.59 


(3P ls)l/2" 


11.52 


(3S ls)l/2+ 


11.62 



TABLE IX. The mass spectrum of Jl° 6c baryons. 
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FIG. 1. The picture of doubly heavy baryon QQq with the colored fields forming the strings between 
the heavy and light quarks, that destroys the pair interactions and involves the additional 'centre-of-mass' 
point close to the centre of mass for the heavy-heavy system. 
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FIG. 2. The spectrum of baryons, containing two 6-quarks: E bb and E bb , with account for the 
spin-dependent splittings of low-lying excitations. The masses are given in GeV. 
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FIG. 3. The spectrum of H+ + and H+ baryons. The masses are given in GeV. 
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FIG. 4. The spectrum of Sj, and baryons without the splittings of higher excitations. The masses 
are given in GeV. 
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FIG. 5. The spectrum of Q® bc baryons with account for the spin-dependent splittings for the low- lying 
excitations. Masses are given in GeV. 
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